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Abstract 

We discuss general positivity conditions necessary for a definition of a 
relativistic diffusion on the phase space. We show that Lorentz covariant 
random vector fields on the forward cone > lead to a definition of 
a generator of Lorentz covariant diffusions. We discuss in more detail 
diffusions arising from particle dynamics in a random electromagnetic 
field approximating the quantum field at finite temperature. We develop 
statistical mechanics of a gas of difi'using particles. We discuss viscosity 
of such a gas in an expansion of the energy momentum tensor in gradients 
of the fiuid velocity. 

1 Introduction 

An approximation of the dynamics of a non-relativistic particle in an environ- 
ment of other particles or random fields by a diffusion process has a long history 
of success [1][2]. A direct generalization of the definition to the diffusion in con- 
figuration Minkowski space does not work because the Minkowski space does 
not have a non- negative scalar product [3] [4] (see the reviews in [5] [6]). One 
could approach the problem of finding a realization of the diffusion by means of 
random dynamics [2] [7] . Kramers diffusion in the phase space can be obtained 
by means of a random perturbation of Hamiltonian dynamics [7] [8] [9] . A rel- 
ativistic generalization of the Kramers diffusion is uniquely determined by the 
requirement that the diffusion process should stay on the mass-shell. Then, its 
generator must be the Laplace-Beltrami operator on the mass-shell [10] [11]. In 
this paper we relax the mass-shell condition admitting a continuous mass spec- 
trum (see some earlier papers on diffusions in four-dimensional momentum space 
[12][13]). In fact, only stable free particles have a definite mass in relativistic 
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quantum field theory. Composite particles and resonances would have a certain 
mass distribution. When we give up the mass-shell condition then the manifold 
of possible rclativistic diffusions is increasing substantially ( a transport theory 
of quantum particles with a continuous mass spectrum has been initiated in [14], 
see also [15]). Any Lorentz covariant positive definite matrix defines a diffusion. 
However, if the diffusion matrix is to depend only on the momentum p^^ then 
its arbitrariness is restricted to a scalar function of . Such diffusions have no 
equilibrium ( they are analogs of the Brownian motion). On general physical 
grounds the equilibrium depends on the frame of reference [16] [17] [18]. We can 
describe the frame by a time-like unit vector -u;^ . If the diffusion matrix can 
depend on both w and p then the set of diffusions is much larger. We deter- 
mine the diffusions resulting from the requirement that they should have the 
covariant (depending on p and w) Jiittner equilibrium distribution [19]. 

The plan of the paper is the following. In the next section we discuss a posi- 
tivity condition resulting from the assumption that the diffusion is defined by a 
probability measure and is rclativistic invariant. An unexpected relation of the 
diffusion matrix to the energy momentum tensor is pointed out. This relation is 
suggesting methods for a construction of rclativistic diffusions. In this section we 
introduce still another method of finding a positive definite and Lorentz covari- 
ant diffusion matrix: the random dynamics. The square of a Lorentz covariant 
vector field defines a positive definite diffusion matrix. If $ is the probability 
distribution of the diffusion process then it defines a current N^^ and the energy 
momentum tensor of a stream of diffusing particles. In sec. 3 we investigate 
whether these (vector and tensor) currents can be conserved in analogy to the 
non-relativistic Brownian motion (in the momentum space). We determine the 
diffusion matrix and the drift of such conservative diffusions. In sec. 4 we briefly 
discuss a relation of the expectation value of the square of a random vector field 
to the random dynamics ( Kubo's diffusion approximation [20] [21]). In sec. 5 
we consider a particle motion in a random scalar and electromagnetic fields. 
We calculate an expectation vahic of the sqiiarc of the vector field defining the 
random dynamics. We show that the diffiision generator determined by the 
random dynamics depends only on the field two-point correlation function at 
coinciding points. The correlation functions at finite temperature depend on the 
reference frame described by the unite time- like vector w^. Using the vector 
w'^ we can define a drift by means of an expectation value of the electromag- 
netic field Ff^i,. We have introduced such a drift in our earlier papers [22] [23] 
as a friction which brings the diffusion to the Jiittner equilibrium on the basis 
of the detailed balance condition (a generalization of the Ornstein-Uhlenbeck 
process). In sec. 5 we show that the diffusion process with a continuous mass 
spectrum can be decomposed into the processes on the mass-shell introduced 
first by Schay [10] and Dudley[ll]. In sec. 7 we introduce some notions of (non- 
equilibrium) statistical physics (energy, entropy, free energy) for a stream of 
diffusing particles and investigate their time evolution and a relation between 
them. In sees. 7 and 8 we initiate a hydrodynamic description of a stream of 
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relativistic diffusing particles in analogy to the well-known description resulting 
from the Boltzmann equation [24] [25] [26] . We show that the diffusion equation 
leads to hydrodynamic equations similar to the ones of the Boltzmann fluid. As 
a consequence we can apply the standard methods of an expansion in gradients 
[25] [26] for an approximation of these equations by a viscous (Navier-Stokes) 
flow. 



2 Some general posit ivity conditions on the dif- 
fusion tensor 

We begin with a general consideration of requirements imposed on the relativis- 
tic diffusion. Explicitly Lorentz invariant evolution can be formulated in terms 
of the proper time. Diffusion can be considered as classical dynamics in a ran- 
dom field. In such an approach the proper time arises in the description of the 
relativistic diffusion. Assume that there is a transition function Pt{x,p; x',p') in 
the proper time r. Then, the diffusion generator (restricted to the momentum 
variables) is 

A = di^'di.d^ + V'd^ = d^idi^'d^ + 6") - d^hi" - {d^di"')di, (1) 

(we denote the derivatives over pi^ by 3^ and over x** by d^). In eq.(l) we 

use two alternative ways to write the generator of the diffusion which will be 
discussed later on. From the definition of the diffusion matrix 

d"- = lim,^o(T-i(pn^) -P^(0))(p'^(r) -p'^(O))) 

= \imr^^T-^ j- dp'Pr{x,p;x' ,p'){p'>' - p^'){p'- - P-), ^ > 

where (..) denotes a probabilistic expectation value. It follows 

a^a.d''^ = lim,_o(aMa<^^"'(p''M -P''(0))(P"W -P'W)) 
= lim^^o r-i((a(p(r) -p(0)))') > 0. 



(3) 



If the process is Lorentz invariant then the generator (1) must be built from 
Lorentz tensors. It is easy to see that a homogeneous (c?'"' independent of p) 
process does not exist if d^^^ is to be built solely from the Minkowski metric 
^nv _ ^]^^ because yy'"^ is not positive definite. Eqs.(l)-(3) hold 

true if r = xo (the coordinate time) and x — x (only the space coordinates 
diffuse) but the Lorentz invariance is not explicit. A formulation of the Lorentz 
invariance is more complicated in this case. 

We may extend the search for relativistic diffusions to generally covariant 
theories. We note that from equations of general relativity the Einstein tensor 

Gfiv = R/xv — 29nvR = Ofjtv, (4) 
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where i?^^ is the Ricci tensor and gij,i, is a pseudoRiemannian metric, satisfies 
the positivity requirement 

a^a.G"" > 0, 

if the energy momentum 9^i, is positive definite . Then, both the Einstein tensor 
and 9 could be chosen as the diffusion tensors. As an example of 9^1, we could 
take the one of an ideal fluid (discussed further in this paper) or the energy 
form of a map ^iM — )• A4 between Riemannian manifolds [27] 

0^.. = gab{<i>)d^,rd.(p\ (5) 

where M is the Minkowski space-time and A4 is the Riemannian manifold 

equipped with the metric gab- 

A general positive definite matrix T^*^ has the representation [28] 

T^"' = J dkG{k)k''k'' , (6) 

where G{k) > 0. In the model (5) we obtain 

G5{k + k') = {gaMr{k)~^\k')), (7) 

if the expectation value (..) is translation invariant ((^ denotes the Fourier trans- 
form of <j)) . 

If G depends only on k and is Lorentz invariant then the argument that a 
positive definite d^^^ does not exist (as discussed after eq.(3)) still applies. We 
cannot define T^'^ by means of eq.(6)using a Lorentz invariant G{k) function of k 
because for such G{k) the integral (6) would be divergent. We shall assume that 
G(fc, w) is a Lorentz invariant function of k and an additional vector w & V+ 
(without the loss of generality we may assume = 1). In such a case the 
tensor T^'^ must be of the form 

T'^'ix) = a{x)w^{x)w''{x) - a;(a;)«" - w^{x)w-'{x)), (8) 

where 77^'' is the Minkowski metric, a and ui are some Lorentz invariant non- 
negative functions. We can express them by T 

a = j dkG{wkf = w^Ti^^w'', (9) 

a; = i y dkG{{kwf - k^) = ^{w'^T^^w'^ - T^. (10) 
If the support of G is in V+ then a > 3a; as 

a = 3uj + j dkGk^. (11) 
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By means of the tensor T^" we can define some new diffusion tensors , e.g., 

d^a = {-VnaTvp + VnpTva " VvpTa^i + VvaTp.p)v'' P'' ■ (12) 

If G is non-negative and is vanishing for fc^ < then eq.(12) defines an admissible 
diffusion tensor (satisfying the inequahty (3)), as we show in sec. 5. 

There is another way to find a positive diffusion matrix. We can find a 
diffusion generator assuming random dynamics. Let A{s) be a random flow (the 
first order differential operator) depending on the proper time s. We consider 

Yr= j dsA{s) (13) 

Then, we define 

= hmt^o y^T-^ = limx^o r'^i dsds' {A{s)A{s') + A{s')A{s)). 

(14) 

If F = a^^d^ then (y+ denotes an adjoint of Y in L?{dp)) 

{Y"") = -{Y+Y) + {a^'d^a'')d^. 

The first operator on the rhs is negative definite (because Y^Y is non- negative). 
It is of the form — d^'^S^Sj^ + 6^9^. Then, it follows from the theory of partial 
differential operators of the second order ( [29], sec. 12) that the matrix d^" sat- 
isfies the positivity condition (3) (because the operator F+F must be elliptic). 
The operator (14) in general is degenerate, i.e., reduced to a lower dimensional 
(as a rule one dimensional) second order differential operator. In general, it is 
only covariant under Lorentz transformations. If the random field is Lorentz 
covariant and we take a mean value of the square in eq.(14) then we obtain a 
non-degenerate invariant negative definite second order differential operator on 
the phase space. This is a good candidate for a generator of the relativistic dif- 
fusion. The randomness may come from a classical approximation of quantum 
fiuctuations. In the next section we discuss some additional conditions imposed 
on the diffusion equation resulting from its probabilistic interpretation. In the 
subsequent sections we show that the methods (6) and (14) are related to each 
other. Moreover, the diffusion generated by F^ of eq.(14) approximates random 
dynamics. 

3 The relativistic diffusion with a continuous 
mass spectrum 

We define the diffusive evolution of observables (functions of position and mo- 
mentum) by the differential equation in the proper time 

drCj> = p^dU + AcP. (15) 
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An evolution of the probability distribution <& in the proper time r is determined 
by the duality 

*r(0) = j dpdxa{p,x)^{p,x)^r{p,x), (16) 

where tr is a Lorentz invariant non- negative function (distribution). The dy- 
namics of the probability distribution <I> in the coordinate (laboratory) time is 
derived from the requirement that $ does not depend on r. This condition gives 
the transport equation 

p^(9;^$ = A+<i> = d^^df^d^*^ + + Q$, (17) 

where A'^ is the adjoint of ^ (1) in L'^{dpa). The functions B and Q can 
be determined from the definition of the adjoint operator and eq.(l). For the 
probabilistic interpretation we need a conservation of the current 
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j dpap''^. (18) 



A Lorentz invariant distribution a in quantum theory of resonances is of the 
form 

(7(/) = j dm^Sijp- - m^)s{m^)9{po), (19) 

where m has the meaning of a mass and 9 denotes the Heaviside step function. 
If a is independent of x then the requirement d^N^ = leads to the equation 

dadp{d°'Pa)-da{B''a)+Qa = 0. (20) 

(under the assumptions that ad, is vanishing at p = and p = oo, so that the 
boundaries do not contribute to the integration by parts formula). 

We assume first that the coefficients in A'^ depend only on p. Prom eq.(20) 
we obtain p°'dap = (because dacr = (j'2pa and from eq.(20) d"''dpa must 
vanish if eq.(20) is to be satisfied for any (j{p'^)) ■ Hence, 

rf'^'^ = Pf'PCpcP'"'. (21) 

where P^" is the projection operator 

= ^/'^ _ (p'^)-'^p^'pr (22) 

If the diffusion depends only on p then from the Lorentz invariance Cpp = 
-VfipP^l- Prom eq.(21) 

d^"" = -jp'^P'"' (23) 

and 

= Xp", (24) 
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where A and 7 > are functions of . Let us note that the matrix d'^^ (23) 
can be derived from eq.(6) when 

with gab — > rj^^. The matrix d of eq.(23) is positive definite if and only if the 
spectral condition arising from eq.(19) is satisfied. 
The current A'''' is conserved if (from eq.(20)) 

Qa = A\a + 2{\u)'p'^ - 6(7cr)y - 127cr. (25) 

We consider also the energy momentum tensor of a diffusing particle 

61'*'' = j dpa'^p^'pr (26) 

The energy momentum is conserved if 

dc^dpid^^Pp^a) - d^{p"b''a) + Qap" = 0. (27) 
Explicitly, with the representation (23) and (24) 

Qa = 5<tA - 1870- + 2(A(7)y - Q{'ya)'p'^. 

We can see that both equations (25) and (27) (for the current conservation and 
the energy momentum conservation) are satisfied if 

7 = ^ n^a{p^). (28) 

In such a case 

A = -K2a(/)p2pA'-a^a, (29) 

(where a(p^) > is assumed to be a regular function at = vanishing at 
infinity). The diffusion generated by the operator (29) is the analog of the non- 
relativistic (Kramers) Brownian motion which also gives a conservation of the 
current and the stress tensor {T^^ = A/'-' , T'*')(the non-relativistic limit of eq.(17) 
with A of eq.(29) coincides with the Kramers diffusion; the diffusion of Schay 
[10] and Dudley [11] has the same non-relativistic limit but does not preserve 
the energy momentum) . If the current and energy momentum conservations are 
satisfied then the transport equation reads 

pi^d^^ = A+^ = -K'^a{p'^)p'^Pi^''d^d„^ + &K^a{p^)pi^dij,^ + (30) 

where 

Q = 12k2 + 6KV-i(aa)y . (31) 

The diffusion (30) has an immediate generalization to an arbitrary metric g'^'^ 
on a pseudoRiemannian manifold M. For this purpose it is sufficient to re- 
place T]'^" in the definition of P'^" (22) by a general metric gi^" and add the 
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geodesic correction V^^p^p" dn<^ (where F's are the Christoffel symbols) on the 
rhs of eq.(30). In such a case the energy momentum tensor (26) of the diffusing 
particles can be put on the rhs of the Einstein equations (4). If the energy 
momentum is not conserved then a varying cosmological constant is needed in 
order to describe the Riemannian geometry of a diffusing matter [30]. 

The diffusion (30) has no equilibrium. An additional friction drift in 
eq.(30) is needed if the diffusion is to equilibrate. The equilibrium state depends 
on the frame of reference [17] [18]. We denote the velocity of the frame by w'^ 
(where w'^Wf^ = 1). Then, the drift is to depend on p^ and on w^. We consider 
a generalization of eq.(30) of the form 

pM^x^ = = -H?a{p'^)p'^P^''d^d^^ + 6K2a(p2)pMa^$ + a'^d^if^^) + Q$ 

(32) 

We have inserted the drift / in eq.(32) in such a way that the current con- 
servation d^N^ = is automatically preserved. With the drift the energy 
momentum conservation fails 

gxg^.. = _ y" dp^f. (33) 

Eq.(33) describes an exchange of energy with an environment (the heat bath). 
The equilibrium is determined by the requirement 

Aj^E = 0. (34) 

can be written in the second form of cq.(l) 

Aj = -a-^d^{K^aa{p'^)p'^P'"'d^ - 2,K^aa{p^)p^' - /"). (35) 

It follows from eq.(35) that if 

f = K'^aap^P^'^d^ In $b - Sn^aap^", (36) 

then eq.(34) for the equilibrium is satisfied. In particular, for the Jiittner dis- 
tribution [19] 

= exp(-/3u;p - Kp'^) (37) 

(where K is an arbitrary non-negative constant and /3 has the meaning of the 
inverse temperature) we have 

Z** = -pK'^aap^P'^'w^ - Sn^aap^". (38) 

The divergence of the energy momentum is expressed again by the energy mo- 
mentum and the conserved current (18)ifa = CT = l. Then 

d^e^"-" = Sk'^N-' - I3k'^{0''^'w^ - w^O"^). (39) 
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Eq.(39) remains true on the mass-shell [31]. We obtain eq.(36) in a model of the 
dynamics of a relativistic particle in a random electromagnetic field in sec. 5. The 
non-relativistic limit of the diffusion (36) coincides with the Ornstein-Uhlenbeck 

process 

dt^ = -pVx$ + dj{dj + p{pj - Wj))^, 

which has a current conservation but no energy momentum conservation (be- 
cause of friction). 

Instead of adding a friction to eq.(30) we first could look for a general 
solution of the positivity requirement upon the diffusion matrix (21). We do 
not know a general solution of the problem. From a study of a particle motion 
in a random electromagnetic field in sec. 5 we obtain 

C"" = -ap^T]"" + b{'n'"'{pwf + p^w^'uf) (40) 

where a and h are non-negative constants. The friction term leading to the 
Jiittner equilibrium (37) will also be derived in sec. 5 from a particle interaction 
with a random electromagnetic field. 

4 Random dynamics 

Explicitly Lorentz invariant relativistic dynamics can be expressed in the proper 
time T [32] 

|^=P», (41) 

-j-=R^[x,p). (42) 
A function W of observables evolves as 

d^W = p^d^W + Ri'd^W = {X + Y)W, (43) 

where 

X = pi'd^. (44) 
The current (18) ($ = T4^) is conserved if 

j dpWd^{R'^(j) = (45) 

Let 

Y{s) = exp(-sX)rexp(sX) = i?^(a; - ps,p)(9^ + sd^), (46) 

where 

Y = R^d^ 

Then, the solution of eq.(43) can be expressed in the form 

Wr = exp(rX)W^/, (47) 
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where 

dsWl = Y{s)W^. (48) 

We assume that i?^ are random variables . In general, we have the cumulant 
expansion for the expectation value 

{W,) = exp(rX) exp ( ds{Y{s)} + \ /J ds ds' {{A{s)A{s') + A{s')A{s))) + . 

(49) 

here 

A{s)=Y{s)-{Y{s)) 

If [Y{s),Y{s')] = and Y is a linear function of Gaussian variables then eq.(49) 
is exact (with no higher order terms). The approach of Kubo [20]- [21] approx- 
imates the random Liouville operator on the rhs of eq.(49) by an expectation 
value of its square. 

The expansion of the dynamics (43) till the second order term reads 

{Wr) =W + Tp'^d^W + Tp^dl /; ds{Y{s))W 

+ /; ds{Y{s))W + i(/; ds{Y{s))fW + i((/; dsA{s)f)W + ... ^^"^ 

In the expansion (50) we find the term of eq.(14). There is also the first 
order differential operator (of the first order in r) which can be defined by 



Ki = lim I ds{Y{s))T- 



Then, there appears in eq.(50) a term which is of the first order in derivatives 
and the second order in time. It can be defined by means of the formula 



K2=\\iTi[ / dsiYis)) - Kit]t-\ (51) 

(F^) -I- K2 determines the diffusion generator (where Y^ is defined in cq.(14)). 
The expansion (49) shows that the diffusion generated by (F^) + K2 is related 
to random dynamics. In fact, Kubo shows [20]that the behaviour in random 
dynamics at times short in comparison to the correlation time goes into the 
difFiisivc T behaviour at times large in comparison to the correlation time. His 
argumcint is equivalent to the rigorous Markov approximation of ref.[9] which 
is using a time rescaling from a microscopic time to the macroscopic time. In 
the next section we calculate the expectation values in eq.(50) for a particle in 
a random electromagnetic field. 

5 Motion in a random scalar and electromag- 
netic fields 

The simplest example of Ri^ in eq.(42) is defined by a scalar field (j) 
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Then, of eq.(14) has the form 



A = T^-'d^d^, (52) 

where T'^" has been calculated in eqs.(6)-(8). We do not develop this diffusion 
any further because in the rest frame {w = (1,0,0,0) and a = 0) it coincides 
with the Brownian motion. We can consider constructed from higher rank 
tensor fields and p^. In this paper we restrict ourselves to 

= TZ^'^p^, (53) 

where TZ^v is the antisymmetric tensor of an electromagnetic field. The electro- 
magnetic current is 

Ju = d^n^^ (54) 

We assume that 

(J") = rw" (55) 
with a certain constant r. We split 

T^f^u = f;," + F^, (56) 

where F"" is the mean value of TZ. F^j^ has zero mean value and is Poincare 
invariant. rWv has the meaning of a constant current. The current (54) results 
from a charge r moving with the velocity w^. We can obtain such a current 
in a finite temperature quantum field theory of interacting electromagnetic and 
complex scalar fields. The density matrix is 

eM-PP"wv) (57) 

where is the four-momentum of the quantum fields and describes the 
moving frame. Then, calculating the expectation value (of the current (54) of 
the quantum scalar complex field (j) interacting with a quantum electromagnetic 
field) we obtain 

= iTr(^eM-PP''^u){^dlct> - ,^9^^)) ^^^^ 
= / rfkfc(^^fc^(exp(/3fcw) — 1)~^ = rw^ 

(with a certain constant r). There is another way to sec that an introduction 
of a current leads to a non-zero expectation value of da-TZ^i, ■ If we add the 
term f Aodx to the Hamiltonian of the quantum electromagnetic field or in a 
covariant and gauge invariant way the term J A^w^dx to the action (here 
is the electromagnetic vector potential) then (7?.^i,,ct) 7^ 0. 

We assume that an average (•) over F is defined which preserves the Lorentz 
symmetry. This means that the two-point function defined by 

{F^,{x)F,p{x')) = G^,,„p{x - x') (59) 
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is a tensor. G/^iy-a-p is symmetric under the exchange of indices x) and 

{ap; x') and antisymmetric under the exchange p, ^ u and u ^ p. We impose 
the Bianchi identities 

^^.,a+^a^,.+^.a,^ = (60) 

on the tensor field 7?.. In terms of the two-point function 

Sae^/^^^G^,;,^ = 0. (61) 
In Fourier transforms eq.(61) reads 

(P^,,{k)F,p{k')) = G^,,„p{k)5{k - k'), (62) 

where G^i,-ap{k) is a tensor which must be constructed from the vectors 

and the fundamental four-dimensional tensors r}fj,p and ep,vpa- Hence,in general 

we could have 

Gf^u;ap{k) = o,i{f]pcrki,kp — t^ppk^k^ + Tji,pkpk^ — T]i,„kij^kp) + do^pivap 
+a2{ripaWukp - ■qppWuka + 'H^pVOpka - rji^aWpkp) (63) 

+a3{1^^laW:,Wp - rinpWi,Wcr + V^pW^Wa - VvaWp.Wp) 

However, the Bianchi identities (61) and the requirement of positivity of the 
probability measure in eq.(59) (sec [33]) lead to oo = 02 = 03 = (in the case 
of a quantum free electromagnetic field at finite temperature we obtain also the 
representation (63) where only oi 7^ 0). We have 



Gp^;ap{x,x') = J dkGp„.cyp{k) cxp{ik{x - x')) 

= J dkGiJz) exp(i/i;(x x ^^{fjpaki/kp i]ppki/k(j + f^i/pk^kfj '^vakpkp^ 
and 



(64) 



{F^^{x)F,p{x'))p = -D^^.,pG{x - x'), (65) 

where 

Dn„;ap = -ripadvdp + nppdj.da - riupdadp. + r]„adpdp (66) 

The two-point function is positive definite if and only if G{k) in eq.(64) satisfies 
the condition 

G{k) > (67) 

and G{k)^0 if P < o[33]. 

It follows from cq.(64) that 

(68) 

where Tp_i, is defined in eq.(6). 
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It is instructive to see the relation between T and the energy-momentum 
tensor O^v for the electromagnetic field. We have 

e.p = - Jry.pF^.F"'^. (69) 

Hence, the expectation value of the energy momentum tensor is expressed as 

(^^.) = 2r^. - \'n^,Tp (70) 

We can show that in eq.(51) 

= {n^.)p''d^ (71) 

and that the rhs of eq.(52) can be expressed as 

( /; ds{Y{s)) - K,t) r-2 = ds ds'd^F^ {x - {s - s')p)p^p'^d^ 

(72) 

Then, taking the limit r — >^ in eq.(72) we obtain 

K2 = {d„n^,)p^p''d^' (73) 

We apply the covariance (65) to calculate the expectation value of the square 
of the Liouville operator. In this way we calculate the term in the expansion 
(49). Then, (A(s) is defined below eq.(49)) 

F2 = lim^^o ^{{Jo A{s)ds)^) = lim.^o ^ /„" ds ds' {A{s)A{s') + A{s')A{s)) 

(74) 

If we apply the formulas (8) then we obtain 
y2 = -2uJ^^,P^"'p^^^ 

+{p + uj){{wpfdf,di' - 2wpwi^pPd^dp + p^w^'wP^^,^p - w^p^d^ - 2wpw''d^) = Ay, 

(75) 

The limit in eq.(74) depends only on DG{0) of eqs.(65)and (68). It does not 
depend on p^. The dependence on of Dij,v-paG{{s — s')p) in eq.(49) is of 
higher than the second order in T(because after integration in eq.(49) over s 
and s' we obtain a factor r^). 

The transport equation determined by the diffusion Y'^ + K2 (eqs.(75) and 
(73)) reads (it is an extension of eqs.(21),(32) and (40)) 

= - aPP-^p^d^ - rp^Pi^^w^ + b{{pwfdi^ - {wp^w^pP + wPpp)dp + p^wf'wPdp)^^ 

(76) 

Here, a = 2uj, b = p + u and the drift (as discussed in sec. 3) resulting from 
eq.(73) is 

fP = -rp^pP'^w^. (77) 

r is an arbitrary real parameter which could be determined from the original 
model of random dynamics(as in eq.(58)). 
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6 Restriction to the mass-shell 

First, we discuss the diffusion at zero temperature corresponding to w = in 
eq. (76). The calculations of (F^) in sec. 5 lead to the following analog of the 
Brownian motion 

dr^ = g^= p^'dl'^ + K^Ap^ (78) 
where is a diffusion constant and 

Ap = -a^P^Va. (79) 

and Q is an operator in the phase space of points (x,p) such that > and 
Po > . The operator — Ap can be defined as a symmetric positive definite 
operator in Li^^dpa) (it is symmetric also in L^{dp)) 

(/i, Ap/2) = D{PdfuP^Pdf2) = J dpa{p)p^P''''dJid,f2. (80) 
Let us consider the direct integral of Hilbert spaces 

H = J dm^s{m^)'Hrn 

where 

(/i J2)m = J d'^pSip^ -m^)h {p)h ip) (81) 
Define the bilinear form 

9{dfi, df2)m = J d^pSip^ - w?)^i^djhdkf2 = (/i, -A™/2)„ (82) 

with 

'yi'' =m'^5^'' +p>p'' (83) 

and 

AS = 7"^5,7^'S*5fc. (84) 

The operator A^ is the Laplace-Beltrami operator on the hyperboloid p^Pfi = 
m?. 7 = det(7jfe) and 7"^ = m^pQ, po = \Jrr? + p^, A; = 1, 2, 3 . The operator 
Aij in Eq.(84) can be expressed in coordinates p as 

A™ = (^i'm2 +p^y)^^ + 3p'^ (85) 

We have 

' dm's{m')g(dh,dh)m = -D{Pdh.v'Pdh)- (86) 



On the Ihs of eq.(86) in each 'Hto we can express the derivatives either by or 
by pj treating po as a function of pj . The relation between the derivatives is 

dj = {djPo)do + dj. (87) 
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The functions fj depend on a four- vector p. On the rhs of eq.(86) we have an 
explicitly Lorentz invariant formula expressed by Lorentz vectors and a measure 
which is a scalar. On the Ihs the Lorentz invariance is not explicit. The Ihs is 
explicitly positive definite whereas the positive definiteness of the rhs is not 
obvious. 

is the generator of the diffusion defined first by Schay[10] and Dudley[ll]. 
The drift (77) can also be restricted to the mass-shell [22]- [23] (it determines 
the exponential speed of the decay to the equilibrium, see [34]). 

7 Thermodynamics of diffusing particles 

We are interested in this section in a thermodynamic description of a system of 
diffusing particles in models which have an equilibrium solution. We find the 
equilibrium solution in the model (76) from the requirement 

{D" - rp'^P'"'w^)<^E = 0. (88) 

We obtain the Jiittner equilibrium (37) if 

r = /3{a-b) = p{ij-a). (89) 

$B is also an equilibrium solution of the model (35)-(36) with r = k^/S. 
Prom eq.(76) we obtain 

gxQtiu = 3^)7^^ + r@i/^w'^ - rG'^'^w^ - 2hw''w^Ni^. (90) 

In the model (35)-(36) we have 6 = a = on the rhs of eq.(90). 

We can introduce thermodynamic notions useful in a description of a stream 
of relativistic diffusing particles. We define the relative entropy current (Kulback- 
Leibler entropy; we follow definitions of our earlier paper [31] concerning the 
model on the mass-shell) 

5^ = Af-i / dpa{p)p''^ ln{N-'^^NE^E^). (91) 

Here, 

N = j dyiNo (92) 

is the charge normalization constant {Ne in the case of the equilibrium). It 
follows that = N^^a^^pQ^^ has the meaning of the probability density (on 
the mass-shell N~^^ is the probability density). It can be shown that 

J dxS°K > (93) 
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and 

d^S^ii = -N-^ j dp<7{p)^a'"'d^ InRd^ InR, (94) 

where R = a^'^ comes from the formula (76) where 

£)M = a^'^d^. (95) 

From positive definitencss of —Aw the momentum dependent coefficients a'"^ 
satisfy the positivity condition a^a^a^" > (see our discussion after eq.(14)). 
It follows that / dxSf^ is a positive function decreasing monotonically to zero 
(it can be interpreted as the entropy of the system plus its heat bath). The 
entropy current of the particle system is 

SI" = -iV-i / dpa{p)p''<!> In (96) 

We have 

d^S" = N-'^ J dpa{p)^ai"'d^ In $9^ In $ - irw^m. 
The entropy is defined as 



j d^S^ 



(97) 
(98) 



Then 



doS = iV-i j dpa{p)^a^"'d^, In ln$ - SrN'^w^, j dxiV". (99) 

We may choose the frame w = (1,0). Then, 

j dxNi, = j dxNo = N>0 (100) 

So, the first term on the rhs of eq.(99) is positive whereas the second is negative. 
We define the free energy 

jr = /3-iAr j dxS% - iV^-^ In(ArAr-i) (101) 

and the energy 

W = y dx(9oo. (102) 

Then, the basic thermodynamic relation (at fixed temperature) 

I3-'^S = W-T (103) 

comes out as an identity. The time evolution of each term in eq.(103) is deter- 
mined by eqs.(39),(91) and (99). From eq.(103) we can see that the decrease of 
OqS in eq.(99) is caused by the decrease of W in eq.(103). 
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8 Fluid velocity and the macroscopic energy mo- 
mentum 

On the basis of the difltusion theory we can develop a hydrodynamic description 
of a gas of diffusing particles. We define the density (a scalar) 

= j dpa{p)^ (104) 

and the mean momentum 

yl^ = pi^y'^N" (105) 

We can write the energy-momentum tensor in the form 

e^'-' = pv^'v" + t"" (106) 

where 

T^''' = j dpa{p){p^' - v^'Xp'' - v-")^ (107) 

Let 

= Waf" (108) 
v^' = nu'' (109) 
Then, u^u'^ = 1. We define the projection operator 

h^^iu) = rj"" - u^vT (110) 

such that u^h^" = 0. In general, we can write 

t'"' = -Uh'"' + a^"" (111) 

(for a gas of free relativistic particles in an equilibrium we have ct'"' = 0). So, 
we write the energy momentum tensor in the form 

61'*'^ == Eu^'u" - Uh^"' + a"" (112) 

where 

E = pn^ (113) 

We shall discuss in more detail the models (35)- (36) and (76) with 6 = (in 
these models fiuid equations have a simple interpretation). As discussed in 
[35] [36] [37] the relativistic fluid equations arc projections of the conservation 
laws. Let us identify the divergence equations (39) with the equations resulting 
from the definition of the energy-momentum tensor (112) (calculating hayd^O^'^ 
in two ways) 
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Let us note that the fluid equations do not depend on the A^'^-term in eq.(90). 
When there is no friction = ) we obtain the same equations as the ones of 
Landau and Lifshitz [36]. The non-relativistic hmit of eq.(114) gives (Euler or 
Navier- Stokes depending on cr'^'^) fluid equations. 

9 Viscosity of relativistic diffusing particles 

The velocity in the transport equation (76) is an arbitrary auxiliary variable. 
In this section we allow to depend on x. We can interpret an equation with 
X- dependent w cither as the diS'usion with friction in a frame moving with the 
local velocity {x) or as the diffusion of a stream of particles in a fluid moving 
with the velocity w^^{x) and playing the role of the heat bath. In the state 

$£; = exp(x(a;) — /3(a;)pii;(x)) (US) 

the particle mean velocity is u'^{x) = 'w'^{x). satisfles the equation A:^^e = 
but the Ihs of eq.(76) is different from zero depending on the derivatives of 
the velocity w, temperature and the density cxp(x(.T)). Then, $£;(x, Xg) 
evolves from an initial time Xq according to the diffusion equation (76). The 
initial and its subsequent time evolution constitute a useful description 
of a state close to the local equilibrium. In this section we discuss the time 
evolution of ^e (115) and the fluid velocity in more detail ((115) is an analog of 
the local equilibrium of the non-relativistic Boltzmann equation discussed in an 
elementary way in [26] and in a more advanced form in [25] ; for the relativistic 
case see [38] [39] and [40]) . 

The mean momentum of the stream in the state is 

PeVe = J dpa{p)pt' exp{-l5pvw''{x)) = UEPEW^ix), (116) 

where 

Pe = j dpa{p)exp{-Pp^w''{x)). 
He is the normalization of v'^ = hew^- We have 

tiePe = f dpa{p)pwex.p{-Ppw) = J dm'^s{m'^)nmPE (H^) 

where rim in 

rimPE = J dpd{p'^ -m^)pwexp{-pp^w''{x)) (118) 

is the value of n calculated on the mass-shell (see e.g. [39]). 

In this section we consider the models (35)-(36) or (76) with 6 = discussed 
in sec. 7. We express the transport equation as the diffusion equation in an 
invariant time variable 

r = x^^w^" (119) 
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( different from r in eqs.(15) and (41); the deterministic proper time satisfies 
T = a;o in the particle rest frame whereas for the time (119) we obtain t = xo 
in the observer's rest frame). We may change the coordinates (a;o,x) (t, x) 
and (pojP) — > {q,p) where 

q = wi'p^ (120) 

In such a case 

p^d^ = qdr + pVx + pf^x^d'^w^dr = qdr+V (121) 

Let us consider 

$ = r$B (122) 

Then, 

pi^d^r = AT- {p^'d^ In ^e)T (123) 

In new coordinates 



qdrV = - Dr - {p^'^^ In $£)r (124) 

The solution of eq.(124) can be expressed by an exponential of q~^{A — V). If 

in the lowest order we neglect the second order space-time derivatives and the 
squares of the first order derivatives of w in F then the solution of eq.(124) reads 

r(a;) = l-q-^ [ dsexp((r - s)q-'^A)p^d^ln^Eis,:x.) (125) 

Taking only the lowest eigenvalue i/ > of —q~^A in eq.(125) we obtain the 
relativistic relaxation time approximation from the formula 

ds exp(-zy(T - s))F{s) ~ u-'^F{t) (126) 

true for large u. Then, the lowest order relativistic relaxation time approxima- 
tion reads 

$-$i3 = i^-iq-V3;j$i3. (127) 

We obtain i' if we can calculate the time evolution (125). We are able to do this 
only in the high energy limit q~^p^ — >■ 0. Such a limit is equivalent to the limit 
m ^ in cq.(76). The time evolution (125) can be calculated exactly [41][31]in 
this limit and the relaxation time approximation is achieved with v = jjn"^. In 
general, in order to obtain a local in time fluid equation we need a local in 
time relativistic approximation to the solution (125). The result of such an 
approximation can be inferred from the relativistic invariance. The relativistic 
invariant (local in time) approximation to eq.(125) must be of the form 

r(p, a;) = 1 + d^w^{x){aop''p'' + aiui'^p'') + a2P^df,x /-, 90-1 

+a3p^d^p + aiw''d^X + a5W^'d^p, ^ 
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where aj{T,pw,p^) are scalars depending on the available Lorentz scalars T,pw 
and p^. In the relaxation time approximation we have 

ao = -{vpwy^ 

and aj = for j > 0. 

Inserting this T in eqs.(18) and (105) we can obtain a formula for the cur- 
rent, from eq.(107) the formula for the energy momentum. Wc could derive the 
general tensor form of these expressions. Subsequently, from the equations of 
conservation of the current and the energy momentum we obtain differential 
equations for the velocity, temperature and the density. The calculations have 
been performed explicitly in the relaxation time approximation to the Boltz- 
mann equation in the classic paper [42]. We would like to concentrate here 
on the appearance of the viscosity (which is of interest for high energy physics 
[43]) in the approximation (125). For this purpose it is sufficient if we restrict 
ourselves to /3 = const and x = const . In such a case 

p($)i'^ = nsPEW^ + cipEw''d^Wfj,{x) + C2PEWf^d^w''{x) (129) 

where Ci and C2 can be calculated from eqs.(18) and (105) as integrals over 
p depending on the functions ao and ai. In the case of the relaxation time 
approximation the constants ci and C2 are determined by the integrals 

s''!--'*" = j drn^sim^) j dpd{p^ - m'^){pw)-^<pEP'"^ ■■■■p'"" (130) 

calculated in [42]. The tensor on the rhs of (130) is expressed by and rj^'^ . 
The energy momentum tensor has the form 

^01^^ dapi^p'^ir - 1)$£; = {E + IlE)wt'{x)w''{x) - rjt^'IlE 
+6o{d^w''{x) + d^w'^ix)) + dirj'^'d^w" - 63wt'{x)w''ix)d^w°'ix) (131) 
-5i{wi^{x)w°'{x)d^w''{x) + w''{x)w°'{x)^^w^'{x)), 

where E and He are the energy and pressure of the Jiittner gas (they are con- 
stants if X and /3 arc constants in cq.(115)). 6j can be calculated if the functions 
ao and ai are known (in the relaxation time approximation we determine 6j from 
integrals of the form (130)). Wc can express u hy w from eq.(129) 



We can see that till the terms first order in gradients we have 



and 



y^fj. — yjij. Cin^^w'^dlw'^. 



w^^ = — cin^^u'^d^u'^. 
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(132) 

(133) 
(134) 



Now, we can replace w in 0'^'^ by u and its derivatives. We have 

ei^" = {E + ILe)u''u'' - rji^^IlE + (j'^j^ (135) 

where 

(136) 

Here 

ai=Si + {E + ILE)cinE^ (137) 

The current conservation and eq.(133) give (when x and /3 are kept constant in 
eq.(115)) 

9;^^" = 0. (138) 
Eq.(114) rewritten in terms of a'^^ (136) with (134) expressed by reads 

eu^a^M^ = -hc,u{u){Sod^h^„dy - (ai - <5o)(w^5p2y-), (139) 

where 

E + nE + /3K^(-3nBCin^i + So - Si). (140) 

In eq.(139) we have omitted the squares of the gradients of u^. This is the 
relativistic Navier-Stokes equation with the shear viscosity So- There is no 
gradient of the pressure and no bulk viscosity in eq.(139) because in the state 
(115) (with constant x and /3 ) in our approximation the fluid is incompressible 
and has constant pressure. 

We can sec that the friction changes the value of + II^;. Apart from this 
minor change we obtain the same hydrodynamics which we could have derived 
from the diffusion theory of sec. 3 with = (preserving the conservation of 
the energy momcntum;in the equation for energy balance the friction would 
appear) . In our definition of as proportional to the current we have chosen 
the Eckart convention [35]. This means that in general the energy momentum 
tensor describes a certain heat flow. We would need ui = So to eliminate; the hait 
flow. We do not expect such an equality to be satisfied for the same reason as in 
Weinberg's discussion [44] of the Thomas model when the general decomposition 
of the energy momentum depends on the definition of the temperature. For a 
complete discussion of the heat flow we would need x-dependent x and j3 in 
eq.(115). If there is no heat flow ( tri = ^o) then in the fluid rest frame we 
obtain the Navier-Stokes term Au^ in the fluid equation (139). We cannot 
formulate a closed hydrodynamic scheme in the sense of Israel and Stewart [45] 
in the approximation which we are applying in this section, in order to check its 
causality and thermodynamic stability [46]. For this purpose we would have to 
discuss the hydrodynamic evolution of the general state (115) to higher orders 
of perturbation. 
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10 Summary 



The assumption of a continuous mass spectrum enables us to work in an ex- 
plicitly covariant way treating all components of the momentum on the same 
footing. We have studied the possible relativistic diffusions cithcir by search- 
ing positive definite diffusion matrices or approximating random dynamics by 
a diffusion. The latter method shows that the relativistic diffusion really arises 
in physical systems. We have studied in more detail the random electromag- 
netic field applied in quantum optics for a description of photons in cavities at 
finite temperature. We identified the non-zero expectation value of the electric 
current as a possible source of friction leading to the Jiittner equilibrium. The 
resulting diffusion may have an application to relativistic streams of particles 
encountered in astrophysics [47] and in high-energy physics [48]. In sees. 7-9 
we have shown how to apply standard methods (widely used in the case of the 
Boltzmann equation) for a statistical and hydrodynamic description of a stream 
of relativistic diffusing particles satisfying the diffusion equation. In this way we 
could relate methods based on the diffusion equation [48] to the ones using the 
hydrodynamic equations in heavy ion physics [37] [49] . The relativistic diffusion 
equation can also supply new ways of approaching the problems of relativistic 
statistical physics and relativistic hydrodynamics by means of methods applying 
differential equations rather than the integral equations of the Boltzmann type. 



References 



[1] 


E.M. Lifshits and L.P. Pitaevskij, Physical Kinetics, Pergamon Press, 1981 


[2] 


S. 


Chandrasekhar, Rev.Mod.Phys.l5,l(1943) 


[3] 


J. 


Lopuszanski, Acta Phys.Pol.l2,87(1953) 


[4] 


R. 


Hakim, J.Math.Phys.9,1805(1968) 


[5] 


C. 


Chevalier and F. Debbasch, AIP Conf.Proc.913,42(2007) 


[6] 


J. 


Dunkel and P. Hanggi,Phys.Rep.471, 1(2009) 


[7] 


N 


Van Kampen, Stochastic Processes in Physics and Chemistry, 3rd Edi- 



tion,North Holland,2007 
[8] Y. Krivolapov and S. Fishman, arXiv: 1203.6895 
[9] H. Kesten and G.C. Papanicolaou, Commun.Math.Phys.78,19(1980) 



[10] G.Schay,PhD thesis, Princeton University, 1961 
[11] R.Dudley, Arkiv for Matematik,6,241(1965) 



22 



[12] M. Rigotti and F.Debbasch, Journ.Math.Phys.46, 103303(2005) 

[13] C. Chevalier and F. Debbasch, J.Math.Phys.49,043303(2008) 

[14] S. Mrowczynski, Ann.Phys.(N.Y.)169,48(1986) 

[15] P. Carruthers and F. Zachariasen, Phys.Rev.D13,950(1976) 

[16] P.J.E. Peebles and D.T. Wilkinson, Phys.Rev.l74,2168(1968) 

[17] R.K. Pathria, Proc.Phys.Soc.88,791(1966) 

J.H. Eberly and A. Kujawski, Phys.Rev.155, 109(1967) 

[18] T. Matsui, B. Svetitsky and L.D. McLerran, Phys.Rev.D34,783(1986) 

[19] F. Jiittner, Ann.Phys.(Leipzig)34,856(1911) 

[20] R. Kubo, J.Math.Phys.4,174(1962) 

[21] R. Kubo, M. Toda and N. Hashitsume, Statistical Physics II. Nonequilib- 
rium Statistical Mechanics, Springer, Berlin, 1985 

[22] Z. Haba, Phys.Rev.E79,021 128(2009) 

[23] Z. Haba, Physica A390,2776(2011) 

[24] C. Cercignani and G.M. Kremer, The Relativistic Boltzmann Equa- 
tion:Theory and Applications, Springer, Berlin,2002 

[25] T.G. Cowling and S. Chapman, The Mathematical Theory of Non-uniform 
Gases, Cambridge, 1970 

[26] K. Huang, Statistical Mechanics, Wiley, New York, 1963 

[27] J. Eells and J.H. Sampson, Amer.Journ.Math.86,109(1964) 

[28] I. M. Gclfand and N.Y.Vilenkin, Generalized Functions, Vol. 4, Academic, 
New York, 1964 

[29] M.E. Taylor, Pseudodifferential Operators, Princeton University Press, 1981 
[30] S. Calogero, J.Cosm.Astro.Phys.ll,016(2011) 
[31] Z.Haba,Mod.Phys.Lett.A25,2683(2010) 

[32] L.D. Landau and E.M. Lifshits, Field Theory, Pergamon Press, New York, 
1981 

[33] Z.Haba, Journ.Phys.A44,335202(2011) 

[34] J.A. Alcantara and S. Calogero, Kinetic and Related Models, 4,401(2011) 



23 



[35] C. Eckart, Phys.Rcv.58,919(1940) 

[36] L.D. Landau and E.M. Lifshits, Fluid Mechanics,Addison-Wesley,1958 
[37] P. Romatschke, arXiv:0902.3663 
[38] W. Israel, J.Math.Phys.4,1163(1963) 

[39] S.R. dc Groot, W.A. van Lccuwcn and Ch.G. van Weert, Relativistic Ki- 
netic Theory, North Holland, 1980 

[40] G.M. Kremer, Continuum Mech.Thermodyn.9,13(1997) 

[41] Z.Haba, arXiv:0911.3126 

[42] J.L. Anderson and H.R. Witting, Physica A74,466(1974), A74,489(1974) 

[43] P. Romatschke and U. Romtschke, Phys.Rev.Lett.99,172301(2007) 

[44] S. Weinberg, Astrophys. J.168,175(1971) 

[45] W.Israel, Ann.Phys.(N.Y.)100,310(1976) 

W. Israel and J.M. Stewart, Ann.Phys.(N.Y.)118,341(1979) 

[46] W.A. Hiscock and L.Lindblom, Phys.Rev.D31, 725(1985), D35,3723(1987) 

[47] J.R. Jokpii,Astrphys.J.146,480(1966) 

[48] R. Rapp and H. Van Hees, in Quark-Gluon Plasma 4, ed. by R.C. Hwa and 
Xin-Nian Wang, World Scientific, 2010; arXiv:0903.10961 

[49] R.Baier, P. Romatschke and A.U. Wiedemann, 
Phys.Rev.C73,064903(2006) 



24 



